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Abstract 

We demonstrate analytically that complex Langevin dynamics can solve 
the sign problem in one-dimensional QCD in the thermodynamic limit. In 
particular, it is shown that the contributions from the complex and highly 
oscillating spectral density of the Dirac operator to the chiral condensate 
are taken into account correctly. We find an infinite number of classical 
fixed points of the Langevin flow in the thermodynamic limit. The correct 
solution originates from a continuum of degenerate distributions in the 
complexified space. 
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1 Introduction 



Random matrix theory and chiral perturbation theory have given deep insight 
into one of the most persisting problems of high energy physics, namely the sign 
problem. The sign problem prohibits direct nonperturbative numerical simula- 
tions of strongly interacting matter with more quarks than antiquarks, since the 
measure on which the Monte Carlo method in lattice QCD seeks to sample the 
most important gauge field configurations is complex valued in the presence of a 
baryon chemical potential (for an excellent review see Ref. [TJ). 

Besides explaining why the quenched approximation fails at nonzero chemical 
potential [2], random matrix theory and chiral perturbation theory have yielded 
highly nontrivial exact results for unquenched QCD at nonzero chemical poten- 
tial n jSlH]. From the viewpoint of the spectrum of the Dirac operator the effect 
of the chemical potential is always dramatic: The chemical potential breaks the 
antihermiticity of the Dirac operator. This drives the eigenvalues off the imag- 
inary axis and makes the fermion determinant complex. In contrast, one of the 
few facts known about QCD at nonzero chemical potential is that the chemical 
potential has very little effect for low temperatures as long as fi is less than the 
scale set by the nucleon mass. This apparent mismatch between the strong effect 
of ji in the fermion determinant and the small effect of /x in the gauge averaged 
fermion determinant (the partition function) has been coined the Silver Blaze 
problem j5]. The exact solutions for the unquenched average density of eigen- 
values from random matrix theory and chiral perturbation theory show how this 
paradox is resolved [6]: as the quark mass enters the band of eigenvalues in the 
complex plane, the average spectral density becomes complex valued and highly 
oscillatory. Since the period of the oscillations is inversely proportional to the 
space-time volume and the amplitude is exponentially large in the volume, the 
oscillations dramatically affect quantities such as the chiral condensate. In fact, 
the oscillations of the eigenvalue density are responsible for the discontinuity of 
the chiral condensate in the chiral limit, in strong contrast to the Banks- Casher 
relation valid for u. = [TJ. 

The standard methods [81418] used to evade the sign problem in QCD (reweight- 
ing, Taylor series, analytic continuation, canonical ensemble, density of states) all 
become extremely hard to handle numerically when the sign problem is severe, 
i.e. when the cancellations due to the complexity of the weight change the parti- 
tion function by a factor which is exponentially large in the volume. In order to 
understand the range of applicability, it has proven useful to study the interplay 
between the sign and Silver Blaze problems. One lesson that has emerged is the 
importance of the phase boundary of phase quenched QCD, i.e. the theory where 
the complex Dirac determinant is replaced by its absolute value, for simulations 
of full QCD at nonzero chemical potential [TU] . 

Complex Langevin dynamics [20l[21] differs from the approaches mentioned 
above in that importance sampling is not used. Instead the field space is complex- 
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ified, which literally opens up new directions to evade the sign problem. Recently 
it has been shown that complex Langevin dynamics can solve the sign and Silver 
Blaze problems in the case of the relativistic Bose g clS , 1 . 6 . , db weakly coupled 
self-interacting complex scalar field at nonzero chemical potential, in four dimen- 
sions [221123] . Even though the sign problem is severe, the phase boundary of the 
corresponding phase quenched theory poses no obstacle. Promising results have 
also been obtained in heavy dense QCD and related models [23] 

Complex Langevin dynamics is not without its problems, though. It has been 
known for a long time that instabilities and runaway solutions can result in a 
lack of convergence of the Langevin trajectories, which necessitates in some cases 
the use of an adaptive stepsize [27J[28] . Recently, it was shown how a straightfor- 
ward implementation of an adaptive stepsize completely eliminates instabilities 
in heavy dense QCD and the three-dimensional XY model at nonzero chemical 
potential [37]. Even when instabilities are under control, it is known that the 
dynamics can convergence to the wrong result [2S]. This problem has recently 
been studied in some detail in the case of simple models with complex noise [38] 
and in the case of the three-dimensional XY model at nonzero chemical potential 
with real noise [39]. Importantly, the conclusion reached in the latter work is that 
the erroneous convergence is independent of the strength of the sign problem. 

In this paper we test the abilities of complex Langevin dynamics against the 
insights obtained about the Dirac spectrum and the sign problem: For QCD in 
one dimension we show that complex Langevin dynamics evaluates correctly the 
contributions from the extreme oscillations of the eigenvalue density of the Dirac 
operator. This provides further evidence that complex Langevin dynamics can 
solve the sign problem, and that the difficulties encountered in Refs. [3EJEH] are 
not intrinsically related to sign problem but have a different origin. This clearly 
distinguishes this approach from the standard ones mentioned above. 

One-dimensional QCD with U(N C ) as gauge group contains many of the fea- 
tures which characterize the sign and Silver Blaze problems in four-dimensional 
QCD at low temperature [30] Q Moreover, it is exactly solvable, which makes 
it an excellent testground for new ideas j3DH33]- In Ref. [30] it was shown that 
the spectrum of the Dirac operator in one-dimensional QCD at nonzero chemical 
potential is located on an ellipse in the complex plane. Since the chiral conden- 
sate can be viewed as the electric field originating from charges located at the 
positions of the eigenvalues, one would naively conclude that the chiral conden- 
sate is zero when the quark mass is inside the ellipse. However, the unquenched 
eigenvalue density on this ellipse is complex and rapidly oscillating: the correct 
chiral condensate, with a discontinuity when the quark mass goes through zero, 
emerges only when all highly oscillatory complex contributions are taken into 

1 Early studies of complex Langevin dynamics can be found in, e.g., Refs. j2M3Q]- Ref. [31] 
contains a further guide to the literature. Other recent work includes Refs. [32ll36| . 

2 On the other hand, for SU(N C ) the sign problem is not severe in one dimension [30] . 
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Figure 1: Sketch of how complex Langevin dynamics evades the sign problem in 
one-dimensional QCD, in the thermodynamic limit n — > oo. The oscillatory lines 
at y — represent the real (full) and imaginary (dashed) parts of the original 
complex weight p(x). The uniform distribution P(x, y) = 1 when p — |/i c | < y < 
fi+\fi c \ and elsewhere represents the distribution sampled by complex Langevin 
dynamics. All y values in this region are equivalent, giving rise to the degeneracy 
mentioned in the title. The original sign problem is severe (mild) when \fj,\ > |/i c | 

Ol < Ia*cD- 

account properly. Even though the Dirac spectrum in four-dimensional QCD is 
spread out into a band, it is exactly the same structure of the oscillations which 
is responsible for chiral symmetry breaking [6]. 

In this paper we establish a first link between the complex oscillations of 
the spectral density of the Dirac operator and complex Langevin dynamics, and 
demonstrate that the sign and Silver Blaze problems in one dimensional QCD 
are solved by complex Langevin dynamics. We encounter a number of surprises 
not seen before. Most importantly, we find that complex Langevin trajectories 
depend on the initial conditions, even in the limit of infinite Langevin time: the 
dynamics is not ergodic. Nevertheless, all trajectories yield equivalent results, re- 
sulting in a degenerate set of stationary distributions in the complexified space. 
In the thermodynamic limit, this continuum of distributions becomes particu- 
larly simple and is sketched in Fig. [TJ Moreover, it is known that classical fixed 
points can play an important role in localizing the dynamics in the complexified 
field space. In most models accessible to an analytical study, only a handful of 
fixed points exists. Instead, in the case studied here, we find an infinite num- 
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ber of stable and unstable classical fixed points in the thermodynamic limit. As 
in Ref. [15] where the factorization method was tested with random matrix re- 
sults, our results demonstrate the usefulness of exact results as benchmarks for 
numerical methods which seek to deal with the sign problem. 

The paper is organized as follows. After a brief review of the relevant ana- 
lytic results in one-dimensional QCD in Sec. [2j we set up and study the complex 
Langevin dynamics problem in Sec. |3j In Sec. S] we show that one stationary 
distribution of the corresponding Fokker-Planck equation can be found analyt- 
ically for all n, but that, surprisingly, this distribution is not realized in the 
actual dynamics. This is explained in terms of classical flow and fixed points in 
Sec. |5j where a larger set of solutions to the Fokker-Planck equation in the ther- 
modynamic limit is given. It is demonstrated analytically that this continuum of 
degenerate distributions yields the correct result for the chiral condensate. Sec. [6] 
contains a brief summary and outlook. 



2 One-dimensional QCD and the sign problem 

We follow closely Refs. [401142] and consider QCD with gauge group U(N C ) in 
one dimension on a lattice with n points. Here n is assumed even throughout 
and is taken to infinity in the thermodynamic limit. The (staggered) fermions 
obey antiperiodic boundary conditions and chemical potential is introduced as 
usual jl6] . We choose the gauge where all link variables are equal to unity, except 
at the final timeslice. The one flavour fermion determinant can then be written 
as 

det[D{U) + m] = det[e n ^ + e~ n ^ + e n »U + e~ ni *U\ (2.1) 

where U is the remaining link variable and \i c is related to the fermion mass m 
via 

m = sinh \i c . (2-2) 

Since there is no Yang-Mills action in one dimension, the partition function has 
the simple form 

dU det[D(U) + m] = s j^t^±M . (2 . 3 ) 
U(N e ) smh(n/i c ) 

One possible way to implement complex Langevin dynamics for this class of 
theories is to evaluate the remaining group integral over the final link variable 
using complex Langevin dynamics. Such an approach has been explored in Ref. 
[21] in the case of SU(3) and excellent agreement with exact results has been 
obtained. Here we wish to make a connection between Langevin dynamics on one 
hand and properties of the Dirac spectrum on the other hand. For that reason we 
take a different route and first cast the chiral condensate as an integral over the 
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eigenvalue density of the Dirac operator. The resulting integral is subsequently 
solved with complex Langevin dynamics. 

The partition function is independent of the chemical potential. On the con- 
trary, the eigenvalues of the Dirac operator D(U), 

X k ,i = 2 (e " +M ~ e n '-^ , (2.4) 

depend on the chemical potential. Here k = l,...,n, and exp(i6i) with I = 
1,...,N C , are the eigenvalues of U. The eigenvalues lie on an ellipse in the 
complex plane, 



sinh(yu) ) ~^ \cosh(/i) ) ^ ^ 

Consequently, the eigenvalue density, 

p(z]fj) = ^ [ dUdet[D(U)+m] Vf(^A M ), (2.6) 

Z JU(N C ) 



'U(N C ) 

'•J~ 

depends on fi. The chiral condensate (normalized with the one-dimensional vol- 
ume), 

E = -— logZ 

n am 

[(1 + N c ) coth [(1 + iV c ) n/ic} - coth (n/i c )} , (2.7) 



cosh(/i 



is, however, independent of [i since the partition function is. Expressing the chiral 
condensate as an integral over the /i-dependent eigenvalue density, 



,2 P( z '^) 



d z ^ (2.8) 

z + m 

the //-independence of E is far from obvious. We note that the condensate E can 
be viewed as the electric field created by the charge density p(z; p). Naively one 
would therefore expect the condensate to be zero when the mass m is inside the 
ellipse in the complex plane, i.e. when \m\ < | sinh(/i)| or equivalently \p c \ < \p\. 
However, this is in contradiction with the known //-independence of the conden- 
sate. This is the Silver Blaze problem in one-dimensional QCD: it illustrates the 
severe problems encountered when lattice QCD is applied for nonzero values of 
the chemical potential. It is thus of practical interest to understand if complex 
Langevin dynamics is able to reproduce this //-independence. 

In one dimension the essential properties of the eigenvalue density do not 
depend on N c . For all N c the eigenvalues lie on an ellipse in the complex plane and 
the eigenvalue density is a rapidly oscillating complex function with a diverging 
amplitude in the thermodynamic limit (when \p\ > \p c \)- We continue therefore 
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Figure 2: The full //-independent chiral condensate (top) in one- dimensional QCD 
along with the contributions from the constant (middle) and oscillating part 
(bottom) of the eigenvalue density. As in four-dimensional QCD, at nonzero \x the 
discontinuity of the chiral condensate at zero quark mass in the thermodynamic 
limit is entirely due to the complex oscillating part of the eigenvalue density 
(compare to Fig. 2 of Ref. [Ej). 



with the case N c = 1, for which the spectral density has a simple analytical 
form [3D]. If the ellipse is parametrized by an angle a, i.e., 

z = - (e ia+ » - e~ ia ~^) , (2.9) 

the eigenvalue density becomes [30], up to an overall constant which cancels in 
Eq. (12~TTT) below, 

p(a;/x) = l- COsh[T ; ( f + ; a)] . (2.10) 
cosh(n/i c ) 

Observe that the eigenvalue density is complex and is highly oscillating when 
| /i | > |/x c |, as illustrated in Fig. [TJ The chiral condensate is then written as 

E= / 7t p0;/^) s 0;/4 (2-11) 

Jo 27T 



with 



/i) = . u - . r~ . u r , (2.12) 
smh(/x + ia) + smh(// c J 
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and is evaluated as 

E = taDh ("^) (2 . 13) 

cosn(/i c ) 

in agreement with Eq. (12. 7p for N c = 1. In the thermodynamic limit the discon- 
tinuity when fi c goes through zero appears and 

lim E = Sg y v (2.14) 

n->oo COsh(/i c ) 

As was first observed in the microscopic limit of four- dimensional QCD [6], 
the strong complex oscillations of the eigenvalue density are responsible for the 
discontinuity of the chiral condensate at zero quark mass in the thermodynamic 
limit. To see this it is advantageous to split the contribution to the chiral con- 
densate into the contribution from the smooth part of the density (the "1"), 



2n da sgn(/i c ) 

£(«;//) = 0(|/x c | - m) — r^^, (2-15) 



2n ' cosh(/i 
and the contribution from the complex oscillating part 



Etf = / ^ [P(«; /-0 - 1] S(a; p) = E - E Q . (2.16) 

This is illustrated in Fig. [2j The smooth part contributes only for large quark 
mass (i.e. \p c \ > \p\) while the oscillating part makes up the chiral condensate 
when |/i c | < In particular, the oscillations of the eigenvalue density are 

responsible for the discontinuity of the chiral condensate when the quark mass 
goes through zero and provide the solution to the Silver Blaze problem. We will 
demonstrate below that complex Langevin correctly evaluates the contribution 
from these oscillations. 



3 Complex Langevin dynamics 

We interpret p(a; /x) as the complex weight, satisfying the usual relation p*(a; /i) = 
p(a; —p*), and write it as 

p{a-p) = \p{a-p)\e w . (3.1) 

The severity of the sign problem can be assessed via the expectation value of 
the phase e l6 with respect to the phase quenched weight. We hence define the 
average phase factor as 

, a, fn™ da p(ot: a) 

(e l6 ) pq = J§ FV ,W , 3.2 

Jo da\p{a;p)\ 
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Figure 3: Expectation value of the phase factor e ie = p(a; p)/\p(a; p)\ with 
respect to the phase quenched weight as a function of p/p c for p c = 1 and 
n = 4,10,16. 

which is shown in Fig. [3] for various values of n. We find the sign problem to be 
mild (absent) when \p\ < \p c \ and severe when \p\ > \p c \. All the dependence 
on p in Fig. [3] emerges from the denominator in Eq. (13. 2\\ , since the numerator 
is /i-independent. The theory with the phase quenched weight therefore has a 
transition at p — p c : our model behaves exactly as QCD in the region where 
< p < ttlb/S, with p c playing the role of m n /2. 

The complex weight p(a; p) is a highly oscillatory complex function, with 
period 2ir /n. From its real and imaginary parts, 

cosh(np) 

Rep(a;u) = l — rcosma), (3.3a) 

cosh(riyU c ) 

Imp(a;/j) = 7 ^\ sin(na), (3.3b) 
cosh(n/i c ) 

we note that the sign problem is severe when the real part of the distribution is not 
positive-definite and, for large n, the amplitude of oscillations grows exponentially 
as exp[n(|/i| — |/i c |)]. On the other hand, when the sign problem is absent, the 
amplitude of oscillations decreases exponentially in the thermodynamic limit. 

We now apply complex Langevin dynamics to study QCD in one dimension. 
We interpret p(a; p) as the distribution and 

S(a;p) = -log p(a;p) (3.4) 

as the complex action. We complexify the angle a — > x + iy. The discretized 
Langevin equations, for general complex noise, read 

x j+ i = Xj + eK x (xj,yj) + eN R r}f, (3.5a) 

Vj+i = Vj + £ K y(xji Vj) + VtNiv), (3.5b) 
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\l = arcsinh m 



1 



u = arcsinh m 



Figure 4: Condensate for n = 4, 10 (left) and n = 16, 22 (right) as a function 
of \x c for fixed /i = 1. The data points are obtained with complex Langevin 
dynamics, while the lines are the exact result (12. 13ft . 



where e is the Langevin stepsize and Langevin time is i9 = je. The drift terms 
are determined by 



K x = -Ref- 
oa 



a—tx+iy 



K T dS 



a^x+iy 



(3.6) 



with the classical drift term 

dS 
da 

The noise satisfies 



sinh[nGu - y + ix)] 

= m w \ TJ—, : — rr- (3.7) 

a^x+iy cosn(n/i c ) — cosh|n(/i — y + tx)\ 



(Vf) = (rj) = (Vfv}>) = 0, (vfvf) = (vWf) = 2 ^"' ( 3 - 8 ) 

with A^r — Ni = 1. Finally, the Fokker-Planck equation underlying this stochastic 
process reads, in the limit that e — > 0, 

d#P(x, y; ■&) = [d x (N R d x - K x ) + d y {Njdy - K y )\ P(x, y; ■&), (3.9) 

where P(x, y; t?) is a distribution in the complexified space, which should be real 
and positive. We specialize to real noise (iVr = 0) in most of the paper, but 
briefly come back to complex noise at the end. 

We have solved the Langevin equations numerically, using a moderate stepsize 
of e = 0.001 and Langevin times up to i9 = 5 x 10 4 . The results for the condensate 
are shown in Fig. H]for four values of n. We observe excellent agreement with the 
exact results, indicated by the lines. Note that the sign problem is severe when 
— 1 < fi c < 1. The discontinuity at fi c — emerges in the thermodynamic limit. 
We have done simulations for various values of fi to confirm that the condensate 
only depends on fi c and not on /i. 
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4 One stationary distribution 

In order to understand why complex Langevin dynamics has no apparent problem 
with the sign and Silver Blaze problems, we first note that fi can be eliminated 
completely from the dynamics by writing y = fi + y, which removes fi from the 
classical drift term. The Silver Blaze problem is therefore trivially solved by the 
complexification. Taking this one step further, we specialize to the case y = fi. 
The corresponding drift term 



is entirely real, such that there is no dynamics in the imaginary direction (in the 
case of real noise) . We can then interpret the remaining Langevin evolution as a 
real process, shifted in the complex plane, with a force K which can be derived 
from an action, 




(4.1) 




(4.2) 



such that K = —c 
can be written a^H 



dS/dx 



It follows that the associated probability distribution 



Ps(x,y) 



Px(x)S(y 



(4.3) 



with 





(4.4) 



where the constant has been fixed by the normalization condition 




(4.5) 



This distribution is related to the original complex weight as 




(4.6) 



and is real and positive. 

This distribution gives the correct chiral condensate, since 




(4.7) 



3 See also Refs. [Tl [38ll48] for similar cases. 
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Figure 5: Expectation value (y) function of (random) initial condition 

^/initial — A 4 f° r fixed /i = 2, /i c = 1 for n = 8, 16 (left) and 32,48 (right). The 
dashed lines indicate (y) = ^initial- 

The solution we found is indeed a stationary solution of the Fokker-Planck equa- 
tion (13. 9 p when Nj = 0, since K y 5(y — fi) = and K = K x {n = 0). To conclude, 
we have found a stationary distribution in the complexified space. If a trajectory 
is initialized at y = fi, we can consider a real Langevin process for which standard 
arguments can be used to demonstrate that the correct stationary distribution is 
reached in the limit of infinite Langevin time. 

Surprisingly, Pg(x,y) is not the unique stationary distribution: trajectories 
initialized with y ^ \x are not in general attracted to y = \x. To illustrate this, 
we have computed the expectation value (y) for a large number of trajectories, 
starting from different initial conditions ^/initial- The resulting (y) is shown in Fig. [5] 
as a function of initial for n = 8, 16 (left) and n = 32, 48 (right), using 100 initial 
conditions randomly distributed around initial = A 4 for each n. Note that we have 
subtracted \x = 2 from both (y) and 7/mitiai- If all trajectories are attracted to 
the stationary distribution Pg, one should obtain (y) = fi independent of ^initial- 
Instead we find that the average value of y is linearly correlated with the initial 
value when ^initial — A* is small and independent of the initial value when ^/initial _ A 4 
is larger and approaches ±/x c . Despite this, all trajectories yield a value for 
the condensate that is consistent with the exact result @ The y value where the 
crossover between the linear dependence on and the independence of the initial 
condition occurs, depends on n. From the numerical results we infer that in the 

4 Taking the average of the 100 initial conditions we find S = 0.6479(2) for n = 8, 0.6483(3) 
for n = 16, 0.6483(4) for n = 32 and 0.6478(3) for n = 48. The exact result is E = 0.648054. 
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thermodynamic limit 



{fi+\tic\ wn en 2/mitiai > // + |// c |, 
2/initiai when [i - |// c | < y init iai < //+ |jU c |, (4-8) 
yu - |// c | when y initia i < pi - \pL c \. 

We conclude that the dynamics is not ergodic. Nevertheless, the expectation 
values of the condensate are consistent within the error with the analytical result 
for all initial conditions. The stationary distribution at y — \i is only realized 

When initial = yU. 



5 Classical flow and degenerate distributions 

In this section we explain the numerical results observed above. We first show 
analytically why the dynamics is not ergodic and subsequently demonstrate that 
a continuum of distributions exist in the thermodynamic limit, all yielding the 
correct condensate. 

The nonergodicity can be understood from the classical flow. We split the 
force explicitly in real and imaginary parts and write 



K - n AE ~ C D K - n A ~° E B (5 V 
(E — AC) 2 + B 2 D 2 y (E — AC) 2 + B 2 D 2 



in terms of 



A = cosh[n(/i — y)}, C = cos(nx), (5.2a) 

B = sinh[?7,(/i — y)], D = sin(nx), E = cosh(n/i c ). (5.2b) 

Classical fixed points are determined by K x = K y = 0. We find 

— n stable fixed points at 

x = {2k + l)n/n, y = fi (k = 0, . . . , n - 1). (5.3) 

— n unstable fixed points at 

x = 2kn/n, y = fi (k = 0, . . . , n — 1). (5.4) 

— 2n points where the flow diverges (K x — 0, K y — oo) at 

x = 2kTr/n, y = fi±fi c (k — 0, . . . , n — 1). (5.5) 
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Figure 6: Flow pattern for \x = 2, \x c = 1 and n = A. The filled/open black dots 
are stable/unstable fixed points, the diamonds indicate a diverging flow. On the 
line y = fi, K y = 0. The arrows are normalized to have the same length. 

The flow patterns are shown in Fig.[6]for n = 4. For larger values of n, the number 
of fixed points and hence the density of regions where the flow changes direction 
increase. The region bounded by ^±/x c is an attractor region: all trajectories will 
end up here, irrespective of the initial y value. The line y = fi, however, is not 
an attractor due to the alternating stable and unstable fixed points, except when 
starting exactly on it. For obvious reasons, we will refer to the region bounded 
by y = H ± A* c as the inside region, while the two regions where \y — (/,] > \fj, c \ are 
referred to as the outside regions. 

The arrows in Fig. [6] are normalized to have the same length. To indicate the 
strength of the force in y direction, we show in Fig. [7] the value of K y at fixed x = 
ir/ '64. We observe that in the thermodynamic limit K y goes to zero in the inside 
region, whereas its magnitude increases linearly with n in the outside region. This 
is confirmed by the following expressions for K y in the thermodynamic limitjf] 

I n sgnfy — /i)e~ n ^ c '"' M_2/ ^ cos(nx) — ^ 0, inside, 
hm K y = < (5.6) 
I n sgn(/i — y) — > ±oo, outside. 



n— >oo 



On the other hand, the force in the x direction goes to zero both on the inside 



5 Exactly on the borderlines y = /i ± /i c , K y diverges for specific values of x, see Eq. (15.51) . 
Away from these points, the force is directed inwards. 
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Figure 7: Force K y in the ^/-direction as a function of y — fi at fixed x = 7r/64 
for various values of n and n — 2, fx c — 1. 



and the outside, 



KKI Ia* y|) sin(nx) — >■ 0, inside, . _ N 

hm K x = ^ _ ; ; (5.7) 

rae rivlM ifi iMcij sin^x) — >■ 0, outside. 



n— too 



This confirms that for large \y\ the flow is attracted to the region \x — |/i c | < 
y < fi + \fi c \ very efficiently. Once inside, the forces vanish exponentially, with 
the rate determined by the vicinity to the boundary at ±/i c . This explains the 
dependence on initial conditions found in Fig. |5j 

In the thermodynamic limit, the forces in the inside region vanish. The 
stochastic evolution then reduces to simple diffusion in a square well bounded 
by y = fi ± fx c (and periodic in x). In the case of real noise the diffusion is 
one-dimensional. In the case of complex noise, the diffusion is two-dimensional. 
It is now straightforward to deduce the stationary solution of the Fokker-Planck 
equation in the thermodynamic limit. We find 

li m P(„) = H when - 1*1 <l/</*+|Mc|, (5 g) 

I elsewhere. 

In the case of real noise, p y (y) depends on the initial conditions as 



"C I ) 



S(y- ix- |/i c | - e) when y initiai >ll+\ll c 
Py(y) = {S(y - yinitiai) when /x - \/x c \ < y mitia \ < /x + \/x c \ : (5.9) 

5{y - ix + \/x c \ + e) when y initial < /x - \/x c \, 
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Figure 8: Normalized condensate S/S exact as a function of n in the case of 
complex noise, N\ > 0. It is seen that in this case complex noise works in the 
thermodynamic limit, n — Y oo. 

where e j 0, see footnote El For complex noise, p y (y) is determined by the 
normalization condition as 

Py(v) = J",- (5-10) 

The resulting distribution is sketched in Fig. [TJ 

This continuous family of distributions all yields the correct value for the 
condensate. To demonstrate this, we write the condensate as 

/dxdy f f dx 

~2 — P(x,y)H(x + iy;ti) = / dyp y (y) / — E(z + iy; fi), (5.11) 

and evaluate the x integral by contour integration. In terms of z — e tx , we find 

dx^, \ f dz -2e-»+v , rio , 

— Z(x + iy;n)= 7T- ( w y (5-12) 

2tt J\ g \=i2m{z-z 1 ){z — Z2) 

with 

z x = e^-^, z 2 = -e^-^y. (5.13) 

In the inside region, \i — \ji c \ < y < fi+ \fi c \, only the pole at z 2 (zi) contributes 
when /i c > (fi c < 0). The result is 

dx , . sgn(/i c ) 

— £(x + u/;/i) = — — — , (5.14) 
Z7T cosn(/x c J 
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for all values of y in the inside region. The remaining integral over y in Eq. (15. lip 
is now trivially performed and yields unity due the normalization condition. We 
conclude therefore that the correct result for the condensate is obtained and that 
the degenerate distributions are all equivalent. 

It is known that complex noise does not work in general (38] and that is also 
what we find here for finite n. However, in this example complex noise can be 
expected to work in the thermodynamic limit, since in that case the dynamics 
takes place in a square well with infinitely-high walls at y — \i ± \i c . This is 
demonstrated in Fig. El where the condensate is shown as a function of n for 
various values of N\. 

6 Conclusion 

We have established a first link between the complex oscillations of the universal 
microscopic spectral density of the Dirac operator and complex Langevin dynam- 
ics. For QCD in one dimension we have shown how complex Langevin dynamics 
correctly evaluates the chiral condensate given the complex and strongly oscil- 
lating unquenched eigenvalue density. The exact solution of the Fokker-Planck 
equation in the thermodynamic limit shows explicitly how the complex Langevin 
method can deal with the severe sign problem present. Surprisingly we did not 
find a unique solution but rather a continuum of degenerate solutions, which 
all yield the correct chiral condensate. This has been shown analytically in the 
thermodynamic limit and demonstrated numerically for finite systems. 

The exact solution presented here offers a direct analytic indication that com- 
plex Langevin dynamics can solve the sign problem. While there exist examples 
where the complex Langevin method is problematic [28 | [38 | l39], we would like 
to stress that the difficulties encountered with complex Langevin dynamics are 
independent of the severity of the sign problem. For instance, in Ref. [39] it 
was demonstrated that the failure of complex Langevin dynamics is caused by 
an apparent incorrect exploration of the complexified field space by the Langevin 
evolution, similar to the case of complex noise considered in Ref. [38]. In the 
cases where the method works well, such as in Refs. [22l|23] and above, the ther- 
modynamic limit poses no obstacle. This is in strong contrast to the standard 
methods (reweighting, Taylor series, imaginary chemical potential and analytic 
continuation), which work well in small volumes but eventually break down due 
the sign problem in large volumes. Our findings therefore strongly encourage 
further studies of complex Langevin dynamics at nonzero chemical potential. 
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